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Abstract 

We study scalar multivariate non-stationary subdivision schemes with a general 
dilation matrix. We characterize the capability of such schemes to reproduce ex- 
ponential polynomials in terms of simple algebraic conditions on their symbols. 
These algebraic conditions provide a useful theoretical tool for checking the re- 
production properties of existing schemes and for constructing new schemes with 
desired reproduction capabilities and other enhanced properties. We illustrate 
our results with several examples. 
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1. Introduction 

In this paper we study multivariate non-stationary subdivision schemes which 
are iterative algorithms based on the repeated application of subdivision oper- 
ators 

S aW -^1{Z S ), f[ fe+1 l =S , aI »]fW, fc>0. (1.1) 

The linear subdivision operators S a ih] 

depend on a dilation matrix M £ X sxs and on the hnite sequences of real 
numbers a^l = {a^, a £ Z s }. If convergent, subdivision schemes are used e.g. 
for designing curves, surfaces, or multivariate functions. 
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The generation properties of subdivision schemes are well-understood and are 
characterized in terms of so-called zero conditions, see e.g. on the 

mask symbols 



These conditions on the symbols determine if the subdivision limit belongs to 
the function space 

EP TiA = sp£m{x~<e Xx : 7 <= V, A e A}, r C %, Ac C s , 

if the starting sequence f ^ in (jl.ip is sampled from a function in this space. If 
for all functions in EP^.a, the subdivision scheme generates exactly the same 
function from which the corresponding starting sequence f ^ is sampled, then 
we say that the scheme reproduces -EPr.A- 

The main goal of our study is to characterize the reproduction property of subdi- 
vision in terms of zero conditions and some additional conditions on subdivision 
symbols. These additional algebraic conditions provide a useful theoretical tool 
that simplifies the analysis of reproduction properties of existing schemes and 
the construction of new schemes with desired reproduction properties. We would 
like to emphasize that the properties of non-stationary subdivision schemes 
could be characterized in the Fourier domain in terms of the properties of the 
associated basic refinable functions, see e.g [H,[l9|, 22, 27, 33, 3J]. Nevertheless, 



we think that it is more advantageous, for construction of new non-stationary 
schemes, to provide such characterizations in terms of simple algebraic condi- 
tions on the coefficients a^, k > 0, of the corresponding refinement equations. 

The reproduction and generation of EP^^a or its subspaces are crucial for model- 
ing objects of different shapes that e.g. are described by polynomial, trigonomet- 
ric or hyperbolic functions. Thus, these reproduction and generation properties 
of subdivision are important in CAGD, motion planning, iso- geometric analysis 
and for studying approximation properties of subdivision schemes }30l |24| . In 
spite of their importance, reproduction properties of non-stationary subdivision 
schemes have not been yet studied as rigorously as it is done in the stationary 
case when the mask a does not depend on the level k of subdivision recursion. 
In the stationary case, characterizations of reproduction of polynomial spaces 
(corresponding to A = {(0, . . . , 0)}) are given in 0, 0, EEl [23] where, in addi- 
tion to zero conditions on the symbol a(z), extra conditions on a(z) and on its 
derivatives at (1, . . . , 1) are established. 

Characterizations of reproduction properties of binary univariate non-stationary 



schemes are given in 12j and we extend them here to the case of multivariate 
non-stationary subdivision with a general integer dilation matrix M. We study 
the so-called non-singular schemes, i.e the ones that generate a zero function if 
and only if f t°l is a zero sequence. The main result of our paper states that a non- 
stationary non-singular subdivision scheme defined by the masks {a[ fc l , k > 0} 
reproduces EPt.a if and only if there exists r € K s such that the symbols aM (2) 
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and their derivatives D~ia^(z) satisfy the zero conditions f|4. 1 5[) and 

s "ft — I 

v^V fe J(v) = |det(M)| ■ v M — H H ((Mr - r)i - j), 

1=1 j=0 

for all v = e~ x ' M ~ (h+1 , -y e r, A e A and k > 0, see Theorem 13771 for details. 

For a convergent non-stationary subdivision scheme, these conditions are suffi- 
cient for reproduction of EPy,k, see subsection 14.2. f I The parameter r € M' s 
above is the so-called shift parameter that determines the parametrization 

t [ £ ] =M~ k (a + T), aeZ s , teR s , k > 0, (1.2) 

associated to a given non-stationary subdivision scheme. This parametrization 
specifies to which grid points t« the newly computed values /„ ' are attached 
at the k— th level of subdivision recursion. The parametrization also influences 
the starting sequences f^. The choice of r does not affect either the generation 
property or convergence of a subdivision scheme, but it affects its reproduction 
properties. In particular, a correctly chosen shift parameter allows to enrich the 
variety of shapes the scheme can reproduce. 

The paper is organized as follows: in Section [2] we introduce notation and re- 
call some known facts about non-stationary subdivision. In Section [31 for a 
non-stationary subdivision scheme, we define the notions of generation and re- 
production of spaces of exponential polynomials and emphasize the difference 
between these notions. Since the univariate case is certainly simpler to follow, 
in Section [4j we first derive the algebraic conditions that characterize the re- 
production of exponential polynomials in the univariate m-ary case, M = m, 
m > 2. Then, in Subsection 14.21 we extend this characterization to the multi- 
variate setting. Finally, in Section [5j we apply the derived algebraic conditions 
and study subdivision schemes of any arity associated to exponential B-splines 
and exponential box splines. We also show the effect of the correct choice of r 
and also the effect of the renormalization of subdivision symbols on the repro- 
duction properties of subdivision. We conclude Section 5 with several univariate 
and bivariate examples. All examples illustrate that our algebraic conditions on 
subdivision symbols make the construction of new schemes with desired repro- 
duction capabilities and with other enhanced properties more efficient. 

2. Notation and Subdivision Background 

We start this section by setting the notation and continue by recalling some 
known facts about non-stationary subdivision schemes. 

- No is the set of natural numbers that includes zero; 

- Multi-indices are denoted by Greek boldface letters a — . . . , a s ) T £ 
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- Vectors are denoted by boldface letters x E R s , or z E C s : 

- In the multi-index notation we have x a — x" 1 • . . . • x" s , at). — a\\ . . . a s \ 
and \a\ = oei + . . . + a s for a E Ng and x E R s ; 

- The product of a real number y E R and a vector x is denoted by y x = 
(yx 1 ,...,yx s ) T ; 

- The scalar product of two s-dimensional (column) vectors x and z is 
denoted by x ■ z = X\Z\ + . . . + x s z s ; 

- The scalar product x T M of an s-dimensional (column) vector x and a 
matrix M is denoted by x ■ M; 

- For vectors x and y, we have e x = (e Xl , • • • , e Xs ) T and (e x ) y — e x y ; 

- Fourier transform / of a compactly supported essentially bounded function 
/ e ioo(R s ) is defined by /(w) = f /(xje-^dx, u E R s ; 

- A sequence of real numbers indexed by If is denoted by boldface letters 
a = {a a , a E Z s }; 

- A sequence of s-dimension vectors indexed by If is denoted by capital 
boldface letters T = {t a , aeZ s }; 

- The space of bounded sequences indexed by If is denoted by loo(Z s ); 

- By D 7 , 7 E Nq, we denote a directional derivative. 

2.1. N on- stationary subdivision schemes 

A non- stationary multivariate subdivision scheme is an iterative algorithm 
with refinement rules 

(w [fcl ) =/r iI =E a L fc WA fcI . « eZS > fc ^°> ( 2A ) 

and it generates the refined data sequence f[ fe+1 l = {fa +1 \ ol E Z s } from 
f [fe] = a £ Z sy_ The dilation ma trix M E Z sxs is assumed to have all 

its eigenvalues greater than 1 in absolute value. Subdivision schemes are based 
on the application of the subdivision operators S a [k] constructed from the so- 
called subdivision mask = {a^, a e Z s } at level k. Each subdivision mask 
is assumed to be of finite support, i.e. supp(a\ k ^) c [0, Nk} s for some positive 
integer Nk- The (non-stationary) subdivision scheme is denoted by S^aM. k>o} 
and is given by 

' Input : fM, {aW, k > 0} 
For k = 0,1,2,..., 
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for some initial data = aeZ s }. 

The symbols of a non-stationary subdivision scheme are given by Laurent poly- 
nomials 

a W(z) = ^ aL fe] 2 a , ze(C\{0}) s , fc>0. (2.2) 

ctgZ a 

Denote to = |det(M)| and by E the set of representatives of Z s /AfZ s . Clearly, 
F contains = (0, . . . , 0) T . Define the set 

5 = { e 2 ™ M ~ T € : ^ i s a CO set representative of Z S /M T 1 S } , (2.3) 

which contains 1 = (1, 1, ... , 1) T . The to sub-masks of the masks are 

{4+M«. « G Z 5 }, eeF, fc>0, (2.4) 
and their symbols (sub-symbols of the masks) are 

4 fe ] (*) = £ ai fc j MQ , £ E E, k > , (2.5) 

respectively. The symbols satisfy 

a [k] (z) = J2 4 fcl («)> fc >°- ( 2 - 6 ) 

For a given parametrization = {ta\ a. E Z s } in (11.21) . a notion of con- 
vergence for S{ a M, fc>o} is established using the sequence {FW, fc > 0} of 
continuous functions F^l that interpolate the data f ^ at the parameter values 
ta\ a £Z S , namely 

F lk] (t^) =/W, Q£Z S , fc>0. (2.7) 

The scheme SV^*] fc>o} applied to initial data f ^ is called convergent, if there 
exists a continuous Zirazi function g { io] (which is nonzero for at least one initial 
nonzero sequence) such that the sequence {F^, k > 0} converges uniformly to 
9m , i-e. 

ff f[ o] = lim 5 aW 5 a[t -n ■••S" a[0 ]f [ ° 1 = lim F^. (2.8) 

The scheme 5/ a [k] ) ^>o} is called weakly convergent, if the sequence {F^, fc > 0} 
in (I2.7|) converges pointwise to g { [ ] € L oc (R s ). 

Due to linearity of each subdivision operator S a [k) , the limit functions g t [o] exist if 
and only if the subdivision scheme applied to the initial data 5 — {5 a o, qgZ s } 
converges to the so-called basic limit function <j> = gg. In this case, 

9m = E /I 01 ^'-")- 

Differently from the stationary case where all masks are the same, in the 
non-stationary setting, one could start the subdivision process with a mask 



5 



at level I > and get a family of subdivision schemes based on the masks 
j a [£+fc] ; > q| ; £ > o. The corresponding subdivision limits are denoted by 
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lim S a ii+k]S a ii+k-i] ■ ■ ■ <S a mf' ', £ > . 



An interesting fact about the compactly supported basic limit functions (f>^ := 
9$\ & > 0, = (j)) is that they are mutually refmable, i.e., they satisfy the 
functional equations 



^2 a$(f> l ' +1] ( M ■ -«). ^ > 0, with a M the £-th level mask. (2.9) 



at 



In this paper, we consider subdivision schemes that are non-singular, i.e. 
they are convergent (or weakly convergent) and, such that g t [o] = if and only 
if f[°l = 0. The following proposition connects the non-singularity of a non- 
stationary scheme with the linear independence of the translates of <f>™, £>0. 

Proposition 2.1. If a subdivision scheme Si a ik] k>o} * s non-singular, then the 
integer translates of <fffl are linearly independent for each I > 0. 

Proof: If the convergent subdivision scheme fe>o} ^ s non-singular, then 

for any starting sequence d = {d a , a G Z s } we have 

lim S alk+e] . . . S aW d = V d a <f>W (■ - a) = 

k— too c — ' 

if and only if d is the zero sequence. I 

Wc formulate the converse of Proposition 12.11 separately as it requires addi- 
tional assumptions on the functions <jffl . 

Proposition 2.2. Let cjffl £ ^(t 8 ), £ > 0, be compactly supported solutions 
of the refinement equations (|2.9[) and such that their integer shifts are linearly 
independent for each I. If there exists L > such that <fffl additionally satisfy 

C||f||oo < II E - f E ^( Z ' S )' £ ^ L ' 

then the scheme Sr^h] k>0) * s non-singular. 



Proof: By 15, Theorem 13], the assumptions on <ft e > guarantee that the associ- 
ated subdivision scheme S'j a [fc] fe>0 i is convergent. The linear independence of 
yields the claim. ■ 
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3. Exponential polynomials and non-stationary subdivision schemes 

In this section, for a non-stationary subdivision scheme, we define the notions 
of exponential polynomial generation and exponential polynomial reproduction. 
We start by defining the space of exponential polynomials on R s . 

Definition 3.1. For T C Ng and A C C s we define 

EP r ,A = span{x~< e X x , 7 £ T, A £ A} . 

The following two observations motivate our interest in the function space 
EP t ,a. 

Remark 3.2. (i) An exponential polynomial p(x) — x~ 1 e Xx E EPp y \ is a 
polynomial, if X — 0, or is an exponential function, if 7 = 0. 7/7 = 
and X E iM. s , then p is a trigonometric function, or a hyperbolic function, if 
7 = and A £ I s . The reproduction and generation of i?r,A or its subspaces 
by non- stationary subdivision are important in CAGD, motion planning or iso- 
geometric analysis. 

(ii) We would like to emphasize that the definitions and the proofs of the 
results of this Section and Section [7] still apply if one works with a subspace 
EPq of EPr t A; where Q C T x A consists of pairs (7, A) for some 7 £ T and 
X £ A. In this case, the algebraic conditions in (|4.2[) . (|4.13l) . (|4. 1 5[) . (|4.21[) and 
(|4.23[) should be checked for corresponding pairs 

(v,7), v=( ei e~^- <fc ^ eE~, (A, 7) £ Q. 

Since most of the properties of a subdivision scheme, e.g. its convergence, 
smoothness or its support size, do not depend on the choice of T[ fc ] = {t„' , a £ 
I s }, these are usually set to 

tW = Ar k a, a. Elf, k>0. (3.1) 

We refer to the choice in (13. ip as standard parametrization. On the contrary, 
the capability of subdivision to reproduce exponential polynomials does depend 
on the choice of t„ and the standard parametrization is not always the optimal 
one. We show in this section that the choice 

t [k] = M -fc( a + T ) ; a£ ^ k > q 5 (3.2) 

with a suitable r £ W turns out to be more advantageous. 
We call the sequence {T^, k > 0} with = {M- k (a + r), a £ Z s } the 
parametrization associated with a subdivision scheme and r £ M s the corre- 
sponding shift parameter. 

Definition 3.3. A convergent subdivision scheme S^k] k>o} * s sa ^ t° be 
SPr^-generating, if there exists a parametrization {T^, k > 0} with T^l = 
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{ta = M k (a + t), a G Z 5 '} and r G R s suc/i i/iai /or every initial sequence 
ft°] = {p(t£ ] ), a G Z s } ; p G PP r . 

lim 5 al4+fc )5 a[t +*-n ...5 aM f [0] G EP rA W > 0. 

We continue by defining the notion of EPr ^-reproduction. 

Definition 3.4. A convergent subdivision scheme S'/ a [*] ) /c>o} * s sa *^ ^° ^ e 
.EPr.A-reproducing, if there exists a parametrization {T^, k > 0} = 
{ta = M~ k {a. + t), a. G Z s } and r G M s suc/i i/iai /or every initial sequence 
fM = {p(tS ] ), a G Z s } ; p G £Pr,A, we fcaue 

lim S a [<+it]5 a [<+*-i] • • • S^mf' ' = p W > 0. 

Note that the generation and the reproduction properties are independent of 
the starting level £ of refinement. 

We define next the step- wise reproduction property of subdivision which is easier 
to check than its -EPr^A-reproduction. 

Definition 3.5. A convergent subdivision scheme Sr^h] fe>o} * s sa id t° be step- 
wise -EPr,A-re-producing, if there exists a parametrization {T^, k > 0} with 
TW = {t^' 1 = M- k (a + t), a G Z s } andrel 8 sucft i/ioi /or </ie sequences 
i [k] = {p(t [ a ] ), a G Z s }, k > 0, andpe EP r , A , we have 

ft fe+1 l =5 aW fW. (3.3) 
It has been already observed in 0] in the stationary multivariate case and 



12j in the binary non-stationary univariate setting that for non-singular 



schemes the concepts of reproduction and step- wise reproduction are equivalent. 
Since these results extend easily to the multivariate case we state the following 
proposition without a proof. 

Proposition 3.6. A non-singular non- stationary subdivision scheme is step- 
wise PP^A-reproducing if and only if it is EPr ^-reproducing. 



4. Algebraic conditions for generation and reproduction of exponen- 
tial polynomials 

In this section we derive algebraic conditions that guarantee PPr,A — reproduction 
by non-stationary subdivision schemes. These conditions arc given in terms of 
the symbols {oJ fe l , k > 0} which are evaluated at the elements of the sets 

V k = {( Vl ,...,v s ) T : v i = e J e-< A - M ~ C * +l) >' ) A G A, e G E }, k > 0. 
For k > 0, we also define the following sets of vectors 

Vl = {(v u ...,v s f : v s =e j e-^ M - lk+1 ^, A G A, £ £-\{l}}. 

To simplify the presentation of the results presented in this section, in sub- 
section [4T] we consider the univariate case first and then extend our results to 
the multivariate setting in subsection 14.21 
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4-1. Univariate case 

In the univariate case, M = m > 2, E = {0, . . . , m— 1} and H = { e 27r ™ ±E : 
e G £J} consists of the m-th roots of unity. From (|2.ip . we get that the m-ary 
subdivision scheme is given by the repeated application of m different rules 

(5 aW fW) =/^+ ] E = E a S+A» eei? ' fc ^°- ( 4A ) 

The structure of the sets Vk, V£, k > 0, is a lot simpler in the univariate case, 
namely, 

V k ={ee- Xm ~ l " +1) : A G A, e G 3}, 

and 

^ = { ee - A "^ (fe+1) : A G A, £G3\{1}}, 
for k > 0, respectively. 

Remark 4.1. We remark that since our goal is to make the multivariate ex- 
tension of univariate results as straightforward as possible, in this subsection 
we use notation less common for the univariate setting. For example, we write 

m—X 

c e , E = {0, . . . , m — 1} instead of Cj. 

4-1.1. Generation of exponential polynomials 

The following result characterizes the -EP^A-generation of a non-singular 
scheme in terms of the so-called zero conditions (|4.2I) . The proof of Proposition 
14.21 in the case m = 2 is given in [3(| Theorem 1] . We give the generalization of 
this result to the case m > 2 in the notation familiar to subdivision audience. 

Proposition 4.2. A non-singular non- stationary subdivision scheme defined by 
the symbols {a^(z), k > 0} is EPr ^-generating if and only if 

D^a [k] (v) =0, 7 G T, v G V£, k> 0, (4.2) 

for VI = { ee - Aro ~ <fc+I) : A G A, e G S\ {1}}. 

Proof. Let I > and A G A. We multiply both sides of the non-stationary 
refinement equation (12.91) by e~ Xm x , and get 

e-Am-^M^) = ^e~ Xm ~ {e+1)a 4 l+1 \mx ~ a ) e -A™-" +1, (™x-«) _ (4 3) 

Set $ e (x) = er Xm ' lx (l)^{x). Then, equation g3J becomes 

= ^ ] e- Xm ~ (l+1)a $e+i(mx - a) , (4.4) 
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or, equivalent ly, on the Fourier side, 

where a [i] (e _i ") = ^ ^e~ wa , w£l The integer shifts of <f>W are EP TA - 
generating, i.e. 

EPy.k C span{^ M ( a; — a) : a £ Z} , 

if and only if the functions <J>£ for all A € A satisfy the so-called Strang-Fix 
conditions 

D^ t (u))\ u=2 ^ = for /3e£\{0}, 7 er, 

and $£(0) ^ 0, see e.g. (35| . 

Therefore, evaluation of (|4.5p at 27T/3, /3e£\{0}, yields 

= $/(2tt0) = m _1 oM fe- Am " tl+1, - 2i,rm " 1 ^ S^-npTim -1 ^. (4.6) 

By Proposition 12.11 the system of functions {$^(- — a) : a £ Z} is linearly 
independent for each £ > 0. By [32|], its linear independence is equivalent to the 
fact that the set {^eC : $ e {£ + ^(3) = for all /3 e Z} is empty. Therefore, 
with £ = 27rm _1 /3 € C, the identities (|4.6[) are satisfied if and only if 

a [i] (v) = 0, veYl, t>°- 

Taking derivatives of both sides of (|4.5|l 

>)<7 ^ ^ ^ 

and evaluating them at 2irf3, j3 £ E \ {0} yields, by induction on 7, that the 
non-singular scheme is -EPr.A-generating if and only if the identities in (j4.2[) arc 
satisfied. □ 

We would like to emphasize that the generation properties of subdivision schemes 
are well-understood. Our interest lies in better understanding of their repro- 
duction properties. 

4-1-2. Reproduction of exponential polynomials 

In this subsection, in Theorem 14.41 we derive algebraic conditions on the 
mask symbols {o\ k \z), k > 0} that characterize the -EPr,A-reproduction prop- 
erty of the associated non-stationary subdivision scheme. We start by proving 
Proposition 14.31 that constitutes the main part of the proof of Theorem 14.41 

Proposition 4.3. A subdivision scheme S ^[k] fc>o} is step-wise EPr ^-reproducing 
if and only if there exists a shift parameter r£R such that 

& %+e ( TO ^) 7 ' v " lf>+£ = v mT - T e- mT+T+s (tut-t- e) 1 ' 

/3SZ 
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is satisfied for all v G V k = {ee~ Am <fc+1> : A G A, e G S}, fc > 0, 7' < 7, 
7 G r, £ G E and e G S. 

Proof: Let 7 G T, A G A and define 

/W = (m- fe (a + r)) 7 -e Am "''( Q+r \ a G Z, jfc > 0. 

By (|4.1| and by definition 13. 5[ the step- wise reproduction of this sequence is 
equivalent to the existence of r G R such that 

/Lt = E a lfi^ eeE > aeZ > k ^°- ( 4 - 7 ) 



Multiplying both sides of by e" Am (fc+1) (™*+»»-+e) we get 

E a S + e - * + -)) 7 e~^-' k+1> ^ 

/3ez 

= (m-MH + E + T)) 1 ^" 1 "!"- 1 ) (4.8) 

or, equivalcntly, 

E a »i+e ( w_fe ( a " ^ + r )) 7 e(- A "^ (fe+1) ^ 2mm " £ )-( m ^ 

_ e -2™i- 1 s-(m/3+e-mT+r) ^ m -( fc +l) ( mQ , _|_ £ _|_ 7 e (- Jl » HH1) -27rmi- 1 e)(mr-r) 

Due to the fact that 

e 2 7 rim- 1 £.m/3 = j for aU £ g £ ^ ^ £ ^ 

we have that (|4.8|) is equivalent to 

E* + e(-- fe («-/^ + -)) 7 ^ +£ 



/3GZ 

= e - mr+T+£ ( m -(W)( mQ + e + r )) 7 (4.9) 
for w G Vfe. Considering both sides of (|4.9p as polynomials in m~ k a with 
(m" fe (a - /3 + r)) 7 = E(Z) (""'H* + (-- fc «) 7 " ■ 

and 

(m-( fc+1 )(ma + £ + r)) T = E( 7 ) (m-( fe+1 )(e + r))" (m- fe a) 7 ^ , 

tj<7 V " / 
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we get that (|4.9[) is equivalent to 

= e - mT+T+e (™- (fc+1) (e + r)y' w mr - r , !)£F t , 77 < 7. (4.10) 
Next, we expand 

(m-( w )(mr|T-mT + e))"= E ( ^ ) (™~ (fc+1) (^ - mr + e)) 7 ' (to^t)^ 7 ' , 

7'<») 

and, similarly, (m~ fc (— /3 + r))' 7 . Thus, (|4.10|l is satisfied if and only if 

E* +£ (-^-^) 7 '^ +e = 

6 -mr+r+e ^ m -(fc+l)( T + £ _ mr ^ 7 ^r-r^ U g Vfc, 7 ' < 7/ < 7. (4.11) 

Multiplying both sides by (— m fe+1 ) 7 , we get the claim. ■ 

To formulate the main result of this subsection we define 

7-1 

qo{z) = l, q 7 (z) = Y[(z-j), 7 eN , zeC. (4.12) 

3=0 

Theorem 4.4. A non-singular subdivision scheme SVyss] fc>o} reproduces EP^ t \ 
if and only if there exists a shift parameter r £ R such that 

, n , ru, , f to • v rnT ~ 7 ' qArxiT — t), for aZZ u suc/i f/iaf e = 1, , , „. 
v D a} ' (v) = < _ " ' , . (4.13) 

I U, otherwise, 

for all veV k = {ee- Am_<fc+1) : A £ A, e € 3}, k > 0, 7 € V. 



Proof: Due to Proposition 13.61 it suffices to prove this statement for step- 
wise -EPr,A-reproduction. We use Proposition 14.31 and show that the step-wise 
-EPr^A-reproduction, is equivalent to the identities (|4. 13[) . To this purpose we 
consider q~ l {rna + e), for fixed 7 € T and e £ E, as a polynomial in ma and 
write 

R. (4.14) 

jj<7 

From (|2.6p we get 



Wa™(z) = E E a ™U 97 ("»a + e)z ma+ ^ j£T. 
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Let v G Vk- By Proposition l4.3l and by (|4.14p . the step-wise -EPr.A-reproduction 
is equivalent to 

e£Er)<f a£Z 
= v mr-r ^ -mr+r+e ^ ( mr - r - £ )" 

eG-E >)<7 

eG-E 

for some r € K. The claim follows, due to 



eG-E 



m, e = 1, 

0, otherwise. 



^.,2. Multivariate case 

In this subsection we give a closer look at the multivariate case. We will 
not repeat results that can be easily extended from the univariate setting by 
simply replacing Z, lor C with Z s , R s or C s , respectively. This is the case of 
Proposition 14.51 which we give without proof. 

Recall that m = \det(M) \ determines the cardinality of E = {e , . . . , e m _i} and 
of 5 = { e 2 ™ M ~ T € : £ is a coset representative of M- T Z S /Z S }. 

Proposition 4.5. A subdivision scheme defined by the symbols {a^(z), k > 0} 
is EPr, a- generating if and only if 

D^a [fcI (v) =0, v G Vl, 7 G T, fc > , (4.15) 

forVi = {( Vl ,...,v s ) T : ^e^'^MeA, eG~\{l}}. 

4-2.1. Reproduction of multivariate exponential polynomials 

The multivariate extension of Proposition l4.3l can appear trivial after reading 
the proof of Proposition 14.61 We believe that would not be the case, if we 
omitted its proof. There are also several crucial differences between the proofs 
of Theorem 14.71 and its univariate counterpart. 

Proposition 4.6. A subdivision scheme S ^g[k\ fc> } is step-wise EPr ^-reproducing 
if and only if there exists a shift parameter r£l s such that 

E & MP+e (M-'Py' ^ MP+£ = V^e-^ (m-^ +1 \Mt -T- £)) V 

/3GZ= 

is satisfied for all v£ Vt = {(^i, • ■ • , v s ) T : Vj = £je~( x ' M +1 ^ , A G A, e G 
S}, k > 0, 7' < 7, 7 G r, £ G E and e G S. 
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Proof: Let 7 € T, A 6 A and define 

f[k] = ( M -*( Q + T ))T . e A.M- fc (a+r) ) Q g ^ fc > g 

By definition, the step- wise reproduction of sequences sampled from the exponential- 
polynomials in EPr,\ is equivalent to the existence of t e I s such that 

C=E4 + A- a€Z'. (4.16) 

Multiplying both sides of d4~TB]) by e ^-M- (H1 '(Ma+Mr+ £ ) we get 



/3ez s 

(M-fW)(Ma + e + r)Y e -^-^ +1 HM-r-r \ (4.17) 



or, equivalently, 
E a M/3+ £ ( M ^( a - # + T )) 7 e(- A - M " (fc+1, - 27riM " T «)-( M ^+ £ ) 

= e -2*iM-THM0+s-MT + T) f M -(k+l) {Ma + £ + T) V £ 



7 (-A-A/-< fc + 1 »-27riAf- T |)-(Af-r--r) 



Now, by properties of scalar products we have M T £ • M/3 = £ ■ M 1 M(3 
£ • f3€Z, and therefore e -2™M- T €-M/3 = L ThuSj (gj^ ig equiva i ent to 

E a M, +e (M- fe («-/3 + x))%^ 



M -(fe+l)( MQ! + £ + r )] v Mt-T j ( 41 g) 

for v e Vfe. Considering both sides of (|4.18[) as polynomials in M~ k ot with 
(M"*(a - /3 + r)y = E(Z) (^H* + ^ 0* ' ? 

and 

(Af-( fe+1 )(Afa + £ + T)) 7 = ^ C 7 ) (Af-( fe+1 )( £ + T)) r '(Af- fc a) T ^, 

»7<7 ^ 77 ^ 

we get that (|4.18j) is equivalent to 

E a M/3 +£ (M- fe (-/3 + r)) ? 'v^ 

£ -Mx +T+£ f M -() i+ l) (e + T) y v Mr-r ) y £ ^ 7? < 7 . (4.19) 
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Next, we expand 

(M-( fe+1 )(Mr + T-MT + e)Y =J2 ( V > ) (aT^t - Mr + £ )) V (M- fc T)"" T ' 

and, similarly, (M~ k (-(3 + t)) v . Thus, (l4~T9)l is satisfied if and only if 
E a ^ +e (-M-^) V v^ +e = 

e -Mx+T+e ^ M -(fc+i)( r + e _ Af-r)) 7 v Mx - T , v G V k , 7' < rj < 7. 
Multiplying both sides by (— l)' 7 , we get the claim. ■ 
Define 

s It — 1 

g (zi ! ...,%) = l, ^...^^nH^-i), 7 GNg. (4.20) 

Theorem 4.7. A non- singular subdivision scheme S^k] k>o} reproduces EPp^ a 
«/ and onZj/ i/ t/iere exists a shift parameter t£1 s smc/i £/ia£ 

[fe] / % _ ( ni ■ v AlT ~ T q 1 (MT — t), for all v such that e = 1, 
1 0, otherwise, 

(4.21) 

for all v£V k = {(«i, . . . ,?; S ) T : ^ = e,-e-( A ' M <fc+1 '^, A G A, e G S}, k > 0, 

7 g r. 



Proof: Due to Proposition 13.61 it suffices to prove this statement for step-wise 
reproduction of sequences sampled from elements of EPp^A- We use Proposi- 
tion 14.61 and show that the step-wise i?Pr,A-reproduction is equivalent to the 
identities (|4~2Tj) . Consider q^(Ma + e) = q^(M k+1 M~ k a + e), for fixed 7 G T 
and e G E, as a polynomial in M a, a G Z s , and write 

<7 7 (A/ fc+1 Af- fe a + £) = ^ c 7 ^ r) (AT fc a)\ ^ei (4.22) 

Let v G Vfe. By Proposition 14.61 and by (I4.22[) , the step- wise reproduction is 
equivalent to 

v^a^v) = J2 E a Ma+ £ ? 7 (MQ+£)v lfa+E 

- E E E 41 +e (A/-^)" v Ma+£ 



EE^( M " (fc+1) (M 

v Mt - t £ - Mt+t ?7 (Mt - t) E e 



r-r-£)l v tfT - T £- M7+T+s 
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for some t € K s . The claim follows due to 



to, e = 1, 



0, otherwise. 



In the case the scheme Si a [*] k>o} i s convergent, but not non-singular, the 
conditions (|4.21[) in Theorem 14 .71 are only sufficient for the E'Pr.A-reproduction. 
Subdivision schemes generating exponential box splines, whose translates are 
linearly dependent, are important examples of such convergent schemes. See 
[34l Theorem 4.3] for conditions that characterize the linear independence of 
exponential box splines. 

Corollary 4.8. A convergent subdivision scheme Sr a [k] k>o} reproduces EPy\ 
if there exists a shift parameter t£I s such that 

~t jyy [k] f \ / to • •v Mt ^ t q^{Mr — t) , for all v such that e = 1, 
1 0, otherwise, 

(4.23) 

for all v e {(«i, • • ■ , v s ) T : v 3 = e j e-^ x ' M ~ lh+1 ^ , A e A, e € E}, k > 0, 7 e T. 



5. Examples and applications 

5.1. Shift factors and interpolator]) schemes 

The first important application of Theorem 14.71 as in the binary univariate 
non-stationary case or as in the stationary multivariate case, is the analysis of 
the reproduction properties of interpolatory schemes, i.e. the schemes whose 
masks satisfy 

It is also of importance to analyse the effect of the mask shifts on the reproduc- 
tion properties of the corresponding schemes. The following results are easily 
obtained by combining [12, Corollary 2] and 0, Proposition 3.4] or 12, Propo- 
sition 3] and 0, Lemma 3.1], respectively. 

Corollary 5.1. A non-singular interpolatory subdivision scheme Sj a [*] k>o} * s 
EPr ^-reproducing only if t = 0. Moreover, the generation of EPy.k implies 
that the scheme reproduces the same space EP^^a- 

Corollary 5.2. If a non-singular subdivision scheme S^ a (k) k>o} reproduces 
EPy,k, then so does the scheme SVbW, fc>o} defined by the symbols b^ k '(z) = 
■ a (k \z), k>0 and [3 £ Z s . 
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5.2. Subdivision schemes for exponential B-splines and exponential box splines 
In this subsection we study the reproduction properties of the subdivision 
schemes associated with the so-called exponential B-splines and box splines. The 
symbols of these schemes even in the non-stationary case play a role of smooth- 



ing factors [151 ] . They also determine the generation properties of subdivision 
schemes. 

5.2.1. Exponential B-splines 

It is well known that in the binary case a non-stationary subdivision scheme 
generates univariate exponential polynomials 

p(x) = x~<e Xx AeC, 7 eN , (5.1) 

if its symbol contains factors of the type 

(1 + Tk z) with a suitable r k G C. 

We show that in the m— ary case such exponential polynomials p are generated 
by non-stationary subdivision schemes with symbols containing factors of the 
type 

(1 + r fe z + r\ z 1 -\ h r™" 1 z m ~ l ) with a suitable r fe e C. 



The following result is a generalization of a result in [15j. We present it here 

as it also allows for derivation of masks of exponential box splines in the next 
subsection and illustrates the generation property of the corresponding schemes. 

Proposition 5.3. Let a non- stationary m—ary subdivision scheme be given by 

a [k] {z) = (l + r k z + r 2 k z 2 + --- + r^ 1 z" 1 " 1 ), r k = e Am ^ +1) , AeC, k > , 

(5.2) 

then its basic limit function is <fi(x) = e Xx X[o.i)- 

Proof: First we observe that a subdivision scheme based on the masks (|5.2[) 

with support size N k — m, whenever convergent or weakly convergent, generates 
a basic limit function supported on [0, fc ~ 1 ] = [0, 1] , see for example [9j. The 
subdivision rules corresponding to the symbols in (|5.2[) are given by 

ti + «\e = r%f*, eeE. (5.3) 

Starting the subdivision process with 6, from (|5.3I) we get that the value of the 
basic limit function (j) at any m-adic point 

k 



/ m J Ej , Ej G E , 



3=1 

is 



k 

^) = IK-i 
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At every non m-adic point x = m J £j, sj G E 1 , we can define 



due to 1 1 — rk | < oo, which ensures the convergence of this infinite product. 



For the point x [0, 1) we set 4>{x) = 0. To show the continuity of <j) at non 
m-adic points we consider a sequence of points 

oo 

{x e , I E N }, with x l = J2 m ~ j £j + TO ~ j£ i> A e E > 

3=1 j=k e +l 

such that kg goes to infinity as £ goes to infinity. Then 



lim \4>(x l ) — 4>{x)\ = lim 

i— yoo i— ¥oo 



II' II' II 

j=l J=l j=fc<: + l 



0. 



The non-uniqueness of the representations of m-adic points 

/c fe— 1 oo 

x = ^ m _J £j = ^ m _J £j + m~ k (ek — 1) + ^ m~ J (m— 1) 

j=l j'=l j=k+l 

with 

oo 

^ m _J '(m — 1) = m~ fe , 
i=fe+i 

makes the analysis of continuity at these points more involved. Thus, we need 
to consider two types of sequences that converge to a; from the left and from 
the right. This implies that is continuous at m-adic points if and only if 



'o • • • 'fc-2 'fe 



fc+i / 



or, equivalcntly, 



= n ^-i 1 = »-r 1 n ^-i 1 = ^ ^-i = • • • = r o m 

j=k+l j=k+2 

Hence, for r = e A / m , A e C, and = e A / mfc+1 we get 



n 



= e 
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Proposition 15 .31 implies that the corresponding subdivision scheme generates 
the exponential polynomials e Xx . Next, we use the algebraic conditions on the 
subdivision symbols in Theorem 14.41 to show how to normalize the symbols 
appropriately to ensure the reproduction of exponential polynomials e Xx . 

Lemma 5.4. Let T = {0} and A = {A}, A G C. Then a non-singular non- 
stationary scheme given by 

oS k \z) = KW{l + r k z + rlz 2 + ---+r™- 1 z m - 1 ), r k = e Am ~ <fc+1) , fc>0, (5.4) 

reproduces EP t .a = {e Xx } if and only if = r^ mT+T and rel. 

Proof: In this case 

V k = { e 27ri ^ m e - Xm ' ik+1> : eeE}, k>0. 

We only need to check the algebraic condition a} k ' (v) — mv mT ~ T for v = r^ 1 = 
e -Am fc+1 ag res ^. Q £ -j-^ conc iiti ons i n (|4.13[) are trivially satisfied. By 
Theorem 14.41 the scheme reproduces SPr.A if and only if for all k > we have 

K^m = qW^ 1 ) = mr t mT+I . 
All these identities are satisfied for r e R and = r k ~ mT+T . ■ 

In the following remark we list several important properties of the exponen- 
tial B-splines. 

Remark 5.5. (i) The scheme associated with 

a [h Kz)= (E^) (E s ^ £ ) (5-5) 

\seE / \e€E ) 

for r k — e Xm <l+1> and Sk = e M ™ <l+1) ; A,/i£ C, k > 0, is convergent, YlK, Udil. 
It has the basic limit function <fi = <fii * (f>2 with 4>\{x) := e Xx X[o,i) an d 4>2{x) = 
e^ x X[0,i)- The function <pi * 4>2 is C , locally an exponential function on [0,1) 
and [1,2), globally supported on [0,2) and is a linear combination of e Xx and 
e^ x . 

(ii) In general, the n-fold convolution <f> — fi\ * (5% * ■■ ■ * /3 n is C"~ 2 , locally 
an exponential function on [J — 1, J), J = 1, . . . ,n, globally supported on [0, n) 
and is a linear combination of the corresponding exponential functions. Such a 
function cj> is a basic limit function of the non- stationary scheme given by n-fold 
products of the symbols in (|5.2[) . Thus, exponential polynomials are generated 
by non- stationary subdivision schemes with factors in (j5.2[) . 
(Hi) By Theorem 4-3] and the compact support of the masks, the exponential 
B-splines satisfy the assumptions of Provosition Wl^ if no two purely imaginary 
A and \i satisfy m~ i (\ — fi) = 2ni£ for I > L, L > 0. Thus, the corresponding 
non- stationary schemes are non-singular. 
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In the case of several exponential factors the space EPp t A that is reproduced 
by the associated non-stationary subdivision scheme is at most of dimension 2. 
The results of Propositions 15.61 and 15.71 are consistent with the observation in 



24j for the case m = 2. 



Proposition 5.6. Let T = {0} and A = {A} 7 A 6 C. Then the non-singular 
non- stationary scheme defined by 

a [k] (z) = KW(l + r k z + r 2 k z 2 + ---+r™- 1 z m - 1 ) n , r k = e Am_(fc+1) , n > 2, (5.6) 

reproduces at most -EPr.A — {e Xx , xe Xx }, if = m 1 ~ n r k mT+T and r S K. 
Proof: We only need to check the algebraic conditions (v) — mv mT ~ T and 

da [h] (z) mr-T a -1 -Am" (fc+1) 

dz z—v 

n (|47T3|) are t 
if and only if 



v a ,, z ' — mv mT t t at v = r,. = e as the rest of the conditions 

dz z—v K 

in (14.13[) are trivially satisfied. By Theorem 14.41 the scheme reproduces EPr,\ 



K^m n = mr~ mT+T 



— 1 Tv-rfcl n-l m ( m_ l) 

r k L K m nr k m n 1 — 



mr k mT+T T. 

2 

The first two identities imply that = m 1 -"r f ; mT+T with r = " (m 2 ~ 1) . 

Next, to guarantee the reproduction of the exponential polynomial x 2 e Xx , the 
algebraic condition v 2 d % — mv mT ~ n ~t(t — 1) should be satisfied. Or, 

° dz A z—v v ' ' 

equivalently, we get 

r k 2 K^ (n(n - l)r 2 k m n - 2 ^— + nm n - l r^m{m - l)(m - 2) j = mr k mT+T T{T-l), 
that becomes 

1 , . , n9 1 . . , , n(m — 1) / n(m — 1) \ 
-n(n - l)(m - l) 2 + -n(m - l)(m - 2) = 1 - ; f V - ' - l\ , 

or, equivalently, 

3(n - l)(m - 1) + 4(m - 2) = 3n(m - 1) - 6 
which can be only satisfied for m = — 1. H 

We continue with the analysis of another case. 

Proposition 5.7. Le£ T = {0} and A = {A,/i} 7 A,/J £ C, A ^ |i. T/ien i/ie 
non-singular non- stationary scheme given by 



VeG-E / VeGB 



(5.7) 
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reproduces at most EPy,a = {e Xx , e^ x }, if 



\e£E 



ana r = n. 



Proof: We only need to ensure that a}- k '(v) — mv mT ~ T for v € {f/7 1 , s^ 1 } as 
the rest of the conditions in (|4.13l) are trivially satisfied. By Theorem 14.41 the 
scheme reproduces -EPr,A if and only if 



K [k] m n fe^r) 



Therefore, the reproduction of EPp y A is possible if and only if 



or, equivalcntly, 



(t- n)(l-m) _ (r-n)(l— m) 
'fe — *fc ' 



which is the case if and only if either r = n or A = fi. For r = n the conditions 
on the first derivative of (z) at z = r^ 1 are satisfied if and only if A = fi. ■ 



Remark 5.8. Increasing the number of different factors in (|5.7[) to 3 would 
lead to the system of equations 



f \ " / \ « 

\ee-B / 



(r-2n)(l-m) 



(T—2n)(l — m) 
W k 



in four unknowns t, r^ — e , Sk — e F and Wk — e' 

Numerical experiments show that these are satisfied when either A = T) = [i or 
only two of A, n, fi are different. The analysis of other possible schemes is not 
one of our goals. 
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5.3. Exponential box splines 

In this subsection we describe the structure of the symbols of non-stationary 
schemes associated with exponential box splines. We use them in subsection 15. 61 
to define a non-stationary butterfly scheme. 

Let M = nl,n> 2. Then m = \detM\ = n s and E = {0, n - 1} S . 
Proposition 5.9. Let 

a [k \z) = ^vlz e , fc>0, r k =e x - M ~ k ~\ A e C s . (5.8) 
Then the basic limit function is (f>(x) = e x ' x X[o,i) s ■ 

Proof: As the symbols (z) have a tensor-product structure, the proof is a 
straightforward generalization of our univariate result on exponential B-splines. 
We only need to observe that (see e.g. [2(| ), even for general dilation matrices, 
every point x in the support of <p has the representation 

oo 

x = ^M- J e i! e 3 eE. (5.9) 

3=1 

The difficulty of studying the continuity of <fi arises only at the points having at 
least one M— adic component, i.e., 

k k—1 oo 

3=1 3 = 1 3=k+l 

where li is the standard I— th unit vector of R s and 

oo 

Y M~ j (n- 1)1 

3 = fe+l 




M 



-fc-i 



(I-M- 1 ) l (n-l)l 



= M~ k (M-I) 1 (n- 1)1 = M" fc l 



-k- 



Remark 5.10. By [35, (4-1)] or by example 6], the scheme defined by the 
masks in (|5-8[) converges weakly to an exponential box spline e Xx X[o,i) s > where 
X[o,i) s * s a characteristic function of [0, l) s . To get the similar result for cf> with 
sets of directions T = {v\, . . . v s }, spanT = R s ; other than E , one just carries 
out the affine transformation y = [vi . . . v s ]x for x € [0, l) s . See lla . example 
6] for the masks of the smoother exponential box splines. 
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5.4- A C 2 binary, non- stationary, dual J^-point subdivision scheme reproducing 
conies 

In this and the following subsections we show that our algebraic conditions 
in (|4.13[) can be efficiently used for constructing new univariate subdivision 
schemes with desired reproduction properties and with enhanced smoothness. 

We start by deriving the non-stationary counterpart of the binary dual 4- 
point subdivision scheme in ljj. We show that, compared to the binary non- 
stationary interpolatory 4-point scheme in [l[ , our scheme reproduces the same 

.Ut 



space of exponential polynomials, i.e. span{l , x, e , e~ } with A € 
{0}, but it is C 2 instead of C 1 . Furthermore, among all existing non-station ary 
binary schemes that are C 2 and reproduce conic sections, see [1(1 11, 12, HI EH 
[25|, [3l|], our scheme turns out to be the one with the smallest support. 
In order to define our scheme, we consider a sequence of masks of the form 

a M - I... c [k] c [k] r [k] r [k] r [k] r [k] c [k] r [k] • • • 1 k > 
A — \ , u, C 3 2 , C 01 , C 2 27 c l,2i c 2,1j c 0,2j c 3,1' u }i fc - u i 



and make use of Corollary 14.81 to determine c\ j, i = 0, . . . , 3, j = 1,2, and the 
shift parameter r <E R such that the space span{l, x, e Xx , e~ Xx } is reproduced. 
Namely, let A = {0, A, -A}, A e (RU iSL) \ {0}, and T = {0,1}. Then, for 
E = {0, 1} and k > 0, we have 



V k = 



_ X2 -(k+i) 



eeE, A e A I 



and 



-A2 



-(fe+i) 



A } A} = {-e 



A2 -()=+i) 



_ A2 -(fc+D 



!}• 



Thus, by CorollaryUand Remark|33 (ii), for Q = {(0, 0), (1, 0), (0, A), (0, -A)}, 
we need to solve the following linear system of 8 equations in 8 unknowns 



' aW(v) =0, Vv e V£, 

„ =0, for v = — 1, 



e ^ 2 - (fc+1) ,l}, 



tfoM (z) I 

=2r, for« = l. 



We get t 



| and 



J*] 


_ M _ 


6{wW) 2 + 2wW - 1 






L 0,2 


— c 0,l — 


64( W W) 3 (2(u;[ fc ]) 2 - l)( W [ fc l 


+ 1)' 




J*l 


_ _M 


10( w W) 2 + 2wW -3 




3 


c l,2 


— c 1;1 — 


64( W W) 3 (2( W W) 2 - l)(w;[ fe ] + 


D + 


4' 


[fc] 


[fc] 


-2(wW) 2 + 2wW +3 




1 


C 2,2 


— C 2,l — 


64(w;[ fe ]) 3 (2(u;W) 2 - 1)( W W + 


D + 


4' 


J*] 


_ J*] . . 


2( U ,W)2 + 2«;[ fe ] + 1 






c 3,2 


~~ c 3,l — 


64(u;W) 3 (2( W W) 2 - 


+ 1)' 
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with 

w [k] = 1^2-' fc + 1 )A/2 +e -2-( fc + 1 >A/2^ 

In conclusion, the non-stationary dual 4-point subdivision scheme we propose 
is defined by the k-level symbol 

« [fel (*) = mzwmt^m^^ {z + 1)3 ( z2 + (4( ™ [fcl)2 - 2)z + l ) ' 

■ ((2( W [ fe ]) 2 + 2 W [ fe l + l)z 2 - (8(u)[ fe l) 4 + 8{wW) 3 + 2)z + 2(w^) 2 + 2w^ + 1 
We observe that 

lim Cq\ = — , lim Cj°\ — 77^, lim c\ k \ = —^z, lim = — t^t, 

k— >oo ' 128 fc— >oo ' 128 fc— ^00 ' 128 k— >oo ' 128 

i.e. when fc tends to infinity converges to the mask of the above mentioned 
stationary dual 4-point subdivision scheme. More precisely, the result in [lil 
Theorem 8] , implies that our non-stationary scheme is asymptotically equivalent 
(of "order" 2) to the stationary dual 4-point subdivision scheme in 17|. This 
property allows us to conclude that our scheme is indeed C 2 . 

5.5. A C 2 ternary, non- stationary, dual J^-point subdivision scheme reproducing 
conies 

In this subsection we derive a C 2 ternary, non-stationary, dual 4-point sub- 
division scheme reproducing the space span{l, x, e Xx , e~ Xx } with A £ (R U iR) \ 
{01. Compared with the ternary, non-stationary, interpolatory 4-point scheme 
in [2j, our scheme will have an additional feature of reproducing conic sections. 
Compared with the ternary, non-stationary, interpolatory 4-point scheme that 
reproduces the same space of exponential polynomials and is C 1 0, 0|, our 
scheme will be C 2 . 

As in the previous subsection, we start by defining a sequence of masks of the 
form 

a W - /. . . n r [k] r [k] r [k] r [k] r [k] r [k] r [k] r [k] r [k] r [k] r [k] r [k] D ■ ■ \ 
a — \ , u, c 03 , c 2 , c Q1 , c 13 , c 12 , e ll , c 2 3 , c 2 2 , c-2,1) ^3,3; 1*3,21 c 3,l' u ' J • 

and we derive the coefficients c[ fc J, i = 0, . . . ,3, j = 1, 2, 3, and the shift param- 
eter t £ R by requiring that the scheme reproduces span{l, x, e Xx , e~ Xx }. Let 
A = {0, A, -A}, A £ (R U iR) \ {0}, and T = {0, 1}. Then, for E = {0, 1, 2} and 
k > 0, we have 

V k = le — e~ A3 : e £ E, A £ A J 

and 

V£= {e^e- A3_<fc+1) : e £ E \ {0}, A £ A j . 
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The corresponding algebraic conditions of Corollary 14.81 and Remark 13.21 (ii) 
with Q — {(0, 0), (1, 0), (0, A), (0, —A)} lead to the following system of 12 equa- 
tions in 12 unknown 

aW(v) =0, Vv E VI, 

a W(v) = 3v 2 r, Vve { e A3 ~ <fc+1 \e- A3 " (fc+1) ,l}, 



< 



da[ J} z) \z=v = 0, w = e"V, £=1,2, 



z= q , = 6t for v = 1. 



Solving the above system we get t = — i, 



[k] 




1 








C 0,2 — 


8wW(2w[ fe l 


- l) 2 (4(u>[ fc l) 2 


-3)( W [ fe l 


+ 1)' 




J*] - 




1 






1 


L l,2 — 


8iyt fc ] (2ioM - 


1) 2 (4( W W) 2 - 


3) (tow + 


D + 


2' 


[fc] 




1 






1 


C 2,2 — 


8wW(2wW - 


1) 2 (4( W W)2 - 


3) (tut*l + 


D + 


2' 


[fc] 




1 








C 3,2 — 


8w[ fe ](2w[ fc ] 


- 1)2(4( W W)2 


-3)(u>[ fc l 


+ 1) 





and 



„[*] _ r [fc] 

-0,3 — c 3,l 



16(tuW) 4 + 16(tuW) s 



24tuW(4(wW)2 - l)3(-4(wW) 3 - 4(wW) 2 + 3wM + 3)' 
[fc] _ [fcj _ 16(w[ fe I) 4 - 16(wt fe ])2 -4tyW - 1 1 

'8wW(2w [fel - l) 3 (2u;[ fc ] + l) 3 (4(«;[ fe ])2 - 3)(u-W + 1) + 6' 
48(u;[ fc l ) 4 + 16(u>M) 3 - 32 (to^ ) 2 - + 1 5 



C 1.3 — C 2,l 



„[*] _„W _ 
c 2,3 ~~ c l,l — 



8w[*](2w[ fe ] -l) 3 (2wW + l) 3 (4(u;[ fc ])2-3)(u>[' £ ] +1) 6' 
[fcj _ [fcj _ 80(w[ fc ]) 4 + 32(w[ fc ]) 3 - 48(tuW) 2 - 12tuW +3 



c 3,3 ~~ c 0,l 



where 



24w[ fc l(4(u>W) 2 - l) 3 (-4(u;[ fc ]) 3 -4(t(;[ fe ])2 + 3u;[ fe ] +3) 



1 ^ e 3-< fc + 1 »A/2 + e -3-< fc + 1 »A/2^j 



2 

The resulting fc-level symbol is 

a [ fe l(z) = -z~*KW (z 2 + z + 1) 2 (z + 1) 

V + (4( W [ fe l)2 - 2)z 3 + (16(u;[ fe l) 4 - 16(tut fe l) 2 + 3)z 2 + (4( W [ fc l) 2 - 2)z + lj ■ 
(16(wW) 4 + 16( W [ fe l) 3 + 3)z 2 + (-64( W W) 6 - 64( W W) 5 + 32( W W) 4 + 
32(^[ fe l) 3 - 12(tuW) 2 - 12tuW - 6)z + 16( W [ fc ]) 4 + 16(tuW) 3 + s) 

where 

K lk] 



24tut*](2tu[ fc l - l) 3 (2wW + l) 3 (4(wW)2 - 3)(wW + 1)' 
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This non-stationary dual 4-point scheme is "order" 2 asymptotically equivalent 
to the stationary dual 4-point scheme with the mask 

lim a [fe] - { 35 1 55 77 9 385 385 9 77 55 1 35 \ 
~ \ 1296' ~16' ^1296' 432' 16' 432' 432' 16' 432' ~1296' ~16' ^1296 J ' 

and its symbol satisfies 

lim a [k] {z) = -z- 6 -i-(z 2 + z + l) 4 (z + l)(35z 2 -94z + 35). 
fc->oo 129o 

By [IH Section 5.2], this stationary scheme is C 2 . Thus, applying [3, Theorem 
8] we deduce that our non-stationary scheme is also C 2 . 

The rest of this section is devoted to the multivariate case. 

5.6. Non-stationary butterfly scheme 

In this subsection, we derive a non-stationary butterfly scheme using the 
results in 0. Let r k = e A2_(fc+1) 7 Ael 2 U iR 2 and k > 0. The symbols of the 
non-stationary butterfly scheme we construct are combinations of the symbols 

a [k] , , { 1 + r k ,izi\ 3 (l + r ka z 2 \ h ( I + r k ,ir k> 2Ziz 2 \ . 
BKJ ll (z 1 ,z 2 )= 1 I I I I I ,j,h,£eN, 

corresponding to the three-directional exponential box splines. In this case, 
M = 21. 

E = {(Q,0f,(l, Of, (0,lf,(l, If} 

and 

S = {(l,lf ) (-l,lf,(l,-lf,(-l,-lf}. 
It is easy to check that the symbols 

= 4 (l r k ^r ka z lZ2 B [ 2 %{z) - 2r k:1 z 1 B l 1 % 3 (z) - 2 r k , 2 z 2 B [ 3 k \ t3 {z) 

- 2r kA r k .2Ziz 2 B [ ^ 31 (z)j, k > 0, 

satisfy (l4~23)) for r = {7 e Np : < 4} and for 

VI = {(-rj&rj^f, (r-^-r-^f, (-r~\, -r^f } . 

Thus, the scheme 5 , { a [t] j fc>0 j is E Pa, r— generating. The scheme is also interpo- 
latory, since the symbols satisfy 

z- 3 z- 3 (>1 (z l7 z 2 ) + W {-z 1 ,z 2 ) + a® ( Zl , -z 2 ) + a® {-z x , -z 2 )) = 4, k > 0. 

Thus, by Corollary 15.11 it is also EP\ y p— reproducing with r = (0, 0) T . We 
observe that our algebraic conditions allow us to simplify the construction of 
the non-stationary butterfly scheme in [29[ and to improve its reproduction 
properties. The smoothness analysis of this scheme is done using standard 
techniques in [l6| . 
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5.7. N on- stationary schemes obtained by convolution 

Another example is of a non-stationary subdivision scheme with dilation 

matrix M = I \ \ I with m — 4. In this case we have 



2 

E = {(0, Of, (1, Of , (1, If, (2, If} , S = {(1, If, (1, -If, (-1, if, (-1, -i) 

since the coset representatives of Z S /M T Z S are {(0, Of, (1, Of, (1, if, (1, 2f }. 
The results in @, [l4| imply that the stationary scheme associated with the 
symbol 

a{z) = \b\z), b(z) = Y,z E , 
is convergent. Thus, by 15|, the non-stationary scheme given by the symbols 
a M(z) = i (>!(*)) 2 , bW(z) = Y,4* e , k>0, 

SEE 



IS 



also convergent, if satisfy |1 — < oo. For example, we let 



fcez + 

g A M £ or some A G M 2 U iR 2 . The conditions in (14. 15)) are satisfied for 

the associated V fe '. Thus, this non-stationary scheme generates polynomials 
{x~ r e X x : 7 G Ng, |7| < 1}. Furthermore, due to (/(o.o^- 2 ) = 1) we get that 
the condition 

4 = aW(v)=4.v M — , v = (r fc - 1 1) r^) ) 

is satisfied if and only if r = (0, 0) T . Moreover, e.g. for 7 = (1, Of, the 
condition 

g = v (i.0)£,(x,0) o [fc]( v j = 4 . y M — q{m (0) = 4 • w Mt - t 

cannot be satisfied for v = (r^ 1 ,r^ 2 ) an d t = (0,0) T . Thus, we can only con- 
clude that the non-stationary scheme reproduces e X x with this shift parameter 

T. 

We show next how to rescale the symbols o\ k \z) appropriately in order to 
enlarge the space of exponential polynomials that is being reproduced. To do 
so, let a^(z) = (b^(z)) , where K^ k \ k > 0, are to be determined from 
our algebraic conditions. For v = (rz\, fu^), we get the system of equations 

4JfW = v MT " T , Mr - t = (ri + r 2 , r 2 f , 

%K [k] = v Mr-x 9(i o)(Mr _ t) = v Mx-r (ri + ^ 
4JC W = v Mr-x g(01) ( Mr _ T ) = v Mt - t T 2 . 

The unique solution of this system is 

K^ = -^ Mt - t and x = (l,lf. 
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Thus, the non-stationary scheme given by the symbols 

a^(z) = \v- k MT+T (&M(z)) 2 , &W(z) = J2 *U e , k > 0, 

eeE 

generates and reproduces the space {x^e^' 31 : 7 e N§, | — y | < 1} for some 
A G M 2 U iM 2 and r = (1, 1) T . 

5.8. \/3— subdivision 

The results of this paper also generalize the results of Q on polynomial 
reproduction of stationary schemes to the case of a general dilation matrix. 
This example shows how to determine the degree of polynomial reproduction 
of the approximating y/i— subdivision schemes given in [28l . page 21] from the 
corresponding mask symbol a(z) instead of the iterated symbol a(ziz^ 2 , z\z 2 x \ 
a(z), as it is done in [7[. We also show how to use affine combinations of these 
schemes to improve their degree of polynomial reproduction, i.e. corresponding 
to A = {0}. 

The dilation matrix in this case is 

M = ( _l _l Y M 2 = -31, 
and the mask symbol is given by 

a(z) = - {z\Zi + zf 1 z 2 ~ 1 + zf x z\ + z^ 2 z 2 + ziz^ 2 + z^z^ 1 ) + 
jjj {zi 1 + z 2 + ziz^ 1 ) + - (z^ 1 + Zi + z^ 1 z 2 ) . 

The associated subdivision scheme satisfies zero conditions of at most order 2, 
see [Hj]. The result of Corollary ITS] yields r = (0, 0) in (fTSJ) . which implies that 
the corresponding scheme reproduces linear polynomials. Since, the mask sym- 
bol satisfies at most zero conditions of order 2, the associated refinable function 
has approximation order 2, see (30| . Note that, similarly, the corresponding r 
is (0, 0) for all approximating V3— subdivision schemes given in 0, page 21]. 

Take next the symbols cij(z), j = 1, 2, 3, 4 of the four approximating subdivi- 
sion schemes in |28l page 21] that satisfy zero conditions of order 3 and consider 
their affine combination 

4 4 

a ( z ) = a j ■ ■ °i(*)> H a j = L 

j'=i i=i 

The resulting scheme still satisfies Condition Z3 and in addition the rest of the 
conditions in Corollary 14.81 for k — 3, if 

Ai = iA 3 + A 4 + 2, A 2 = -|a 3 - 2A 4 - 1, A 3 ,A 4 eIR. 



28 



Surprisingly, for any such choice of the parameters Ai,...,A4 the resulting 
scheme is given by 

a(z) = l-hzi 2 +Zi 2 z%+z%+zl+zlz2 2 +Z2 2 )+^(zi 1 +Zi 1 z 2 +z 2 +z 1 +z 1 z 2 ; 1 +z 2 ; 1 ) 
9 9 

and defines the interpolatory scheme considered in [H, page 17]. Thus, by (30l |. 
this scheme has approximation order 3. 
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